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The problem of explanation of the critical angular velocity Q, c when the formation of a vortex in 
the stirred Bose-Einstein condensate becomes energetically possible, is considered in the framework 
of the variational approach. The origin of smallness of the calculated fi c in comparison with the 
measured values which takes place for pure quantum state with the unit angular momentum per 
condensed particle, is uncovered. The agreement with the measured fi c is achieved upon admitting 
a small admixture of the zero angular momentum state in the wave function of the one-vortex 
quantum state prepared after stirring. The portion of this admixture amounts to 10—13% of the 
total condensed atoms. Possible test of this hypothesis is proposed. 

PACS numbers: 03.75.-b, 03.75.Fi,05.30.Jp 

Bose-Einstein condensate (BEC) discovered in the trapped clouds of alkali atoms [[[) are proved to be an excellent 
test site of fundamental concepts of quantum physics of the systems consisting of macroscopically large number of 
particles . One of the most intriguing features of such systems is the quantized vortex. Recently, the formation 
of such vortices in the trapped BEC was observed in two different situations. First is the two-component BEC [[|, 
where the vortex state is created via interconversion between two hyperfine states. The second one is realized by the 
stirring of BEC with the toggled laser beam. This toggling beam creates a small axial asymmetry of the trap potential 
which is rotated slowly with the angular velocity f2. It was found that for f2 exceeding a definite critical value f2 c , the 
images of the space distribution of BEC atoms after the ballistic expansion reveals a visible signature of the vortex 
j|, and even the lattice of the vortex array was observed in some situiations The method of Ref. Q directly 
corresponds to the classical experiment with the rotating bucket ||. As was pointed out in Ref. Q, the measured 
critical angular velocity of the formation of a single vortex in the stirred BEC appears to be notably larger than that 
predicted theoretically § . 

Recently, the efforts aimed at explaining larger critical angular velocity were undertaken in Ref. 0. The purpose 
of the present note is to propose another explanation of the critical angular velocity f2 c of the stirring of BEC when 
the formation of the single vortex becomes energetically possible. To this end the energy of different configurations 
of the BEC atoms is calculated. The reason of smallness of the calculated £l c as compared to the measured one is 
revealed. It is shown that the agreement with the measurements can be achieved by assuming the admixture of the 
vortex-free state with zero angular momentum, in the wave function of the final state prepared after the stirring and 
containing the visible vortex. 

As is known all basic properties of BEC in diluted gases of alkali metals are described by Gross-Pitaevskii (GP) 
equation Q which has the form of Schrodinger equation added with the nonlinear term arising due to the short range 
interaction characterized by the single parameter - the scattering length. Since it is the equilibrium energy of the 
BEC gas that is of the main concern here, the GP energy functional 

E = /d 3 x{|l|V^ 2 + yKri+^z 2 )H 2 + ^H 4 }, (1) 

is used instead of GP equation. In the above equation, — x 2 + y 2 , m is the mass of an atom, u±, uj z are, respectively, 
the transverse and longitudinal frequencies of the oscillator-like potential modeling the axially symmetric trap, a is 
the scattering length. Also, ip is the condensate wave function normalized according to the condition 

N = J d 3 a#| 2 , (2) 

N is the number of condensed atoms. The trap parameters N, cu±, and lo z are specified as follows. The first set 
referred below as the set A is |Q 
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N = (1.4 ± 0.5) x 10 5 , — = 219 Hz, — = 11.7 Hz, 
2n 2ir 



(3) 



and the corresponding critical angular velocity of stirring is fl c /2n — 152 Hz M. The second set referred below as the 
set B is H 



N = (3.7 ± 1.1) x 10 5 , 
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and the corresponding critical angular velocity of stirring is £l c /2ir = 115 Hz M. 

Since the accuracy of determination of the number of condensed atoms is about 30 %, and the energy of BEC is 
scaled as E oc N 2 ^ 5 [see Ref. Q and Eqs. (0), flTl| ) below], it is unnecessary to bother to calculate E with the accuracy 
better than 10%. So, one may hope that the variational calculation []l(| of energy instead of full numerical solution of 
GP equation will be sufficient. As will soon become clear, it is convenient to take the trial wave function in the form 
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where < k < 1; R K and z K are variational parameters. Notice that k = 0, 1 corresponds to, respectively, pure 
vortex-free state and the state with the singly quantized vortex placed at the center of the trap, while intermediate 
values of k correspond to arbitrary mixture of the above states. Introducing the dimensionless parameters p K and 
z K according to the relations R K = (fi/mwi) 1 ' 2 p K and z K — (h/mtUz) 1 ^ 2 £ re , one can find from Eqs. (|l]) and (|J) the 
energy per condensed atom: 
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The mean value of angular momentum in the quantum state with the wave function Eq. ([|) is 
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As is known [£ 11 1, the condition of thermodynamic possibility of the vortex formation in the system rotated at the 
angular velocity SI can be formulated as AE — (L z )Cl < 0, where AE is the energy difference between the states with 
vortex and without it, so that the critical angular velocity is defined as f2 c — AE/(L Z ). Let us evaluate fi c in the 
framework of variational approach. 

The values of variational parameters p K and C, K can be found from the condition of the minimum of the BEC energy 
Eq. (p) which is reduced to the following equations: 
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First, let us consider the problem of the critical angular velocity in the approximation when the kinetic energy of both 
the transverse and longitudinal motion can be neglected. This is the Thomas-Fermi (TF) limit The solution of 
Eq. (||) in this limit looks as 
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Direct numerical evaluation shows that at, say, k 
evaluated in this limit for the set of parameters Eq. i 



the parameters pJ F = 2.16 (2.63) and (£ ¥ = 9.36 (10.7) 
[respectively, (0)], coincide within the accuracy 1% with those 
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found from Eq. ([}]). Hereafter, when doing the specific numerical evaluations, we take the scattering length a — 5.77 
nm for 87 Rb atoms fl^| , and the above two sets of BEC parameters from Eqs. (||) and The corresponding energy 
in TF limit is found to be 
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One can see that in TF limit the energies of pure vortex- free state (n = 0) and the state with the single vortex (k = 1) 
are equal, hence tt c found from relation 



Q c (k) _ E(k) - E(0) 
2?r " 2ttTiNk, ' 



(12) 



[see [|| and Eq. (g)], at K = 1, vanishes in TF limit. Thus, one should take the kinetic energy of the BEC cloud into 
account. Here this is done by the numerical solution of Eq. (||) for the two values n = 0, 1. The result is pa = 2.19 
(2.64), Co = 9.28 (10.66), and p x = 1.61 (1.90), (i = 9.06 (10.53), in the case of the experimental conditions Eq. (§), 
([|), respectively. One then finds = 65.9 Hz, and j2n = 36.1 Hz, where upper indices refer to the two above 

sets of the trap parameters. The smallness of these values as compared to the experimentally measured is due to 
smallness of the energy difference between the pure BEC states with k = and 1. 

To reconcile the result of calculations with the measurements in the present approach, one should have in mind 
that, in fact, the density of atoms does not vanish in the central dip Q]. The authors of Ref. [Q propose three possible 
reasons for this: (i) oscillations of the vortex filament, (ii) the presence of non-condensed atoms, and (iii) insufficient 
resolution of the imaging optics as compared to the vortex core radius of the BEC cloud. Here we propose the fourth 
possible reason and admit that the quantum state of the BEC cloud after the stirring is the superposition of pure 
quantum states with the angular momenta per particle L z /N = and h, that is, admitting k ^ 1 in the wave function 
Eq. (Q) of final quantum state. The presence of BEC atoms in the state with zero angular momentum explains in a 
natural way a nonzero density in the central dip. 

Solving Eq. (|^) numerically, one can find the energy dependence on k and calculate the critical angular velocity 
from Eq. (|l2|). The result of this calculation is shown in Fig. [|. Then fitting the calculated critical angular velocity 
Eq. (|l|) to the experimental values is possible if the portion k of the number of BEC particles in the state with the 
unit angular momentum per particle amounts to, respectively, 



= 0.87, k b = 0.90. 



(13) 



The latter values, in view of Eq. (||), give (Lz) /N = 0.87?i, respectively, O.QTi, and do not contradict to the measured 
magnitude Q of the mean angular momentum of the condensate. Eq. (|l^) means that the admixture of the atoms in 
the zero angular momentum state in the case of the trap parameters Eq. (|^) [(^)] amounts to 0.13 (0.1), respectively. 
Notice that despite essential difference in the experimental trap parameters Eqs. (||) [f| and (jij) ||, the portion of the 
number of atoms in the zero momentum state needed to explain very different observed critical angular velocities in 
the final state prepared after the stirring, turns out to be practically same. To be more precise, the 30 % accuracy 
of the determination of the number of atoms in the condensates reported in Refs. QJ^] implies, as is explained earlier 
in this paper, approximately 10% uncertainty of the calculation of energy per condensed atom which is translated to 
approximately the same uncertainty of calculation of k, while the central values of the calculated k differ in the above 
experimental conditions by 3 — 4%, which is well below their estimated uncertainty. 

Could such proposed feature of the BEC wave function as the presence of the portion 1 — k of zero angular momentum 
condensed atoms be tested in experiments? Let us discuss this issue. As is known 00], the presence of the vortex is 
detected through the visualization of the images of BEC cloud obtained after its ballistic expansion. One can obtain 
the spatial distribution of the BEC atoms after this expansion with the usual quantum mechanical method upon 
finding the wave function Eq. (|^) in the momentum space form, propagating it forward in time freely, then finding its 
resulting coordinate space form. The resulting spatial distribution appears to be 
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One can see that the dependence of the spatial distribution of the BEC atoms on the azimuth angle <f> is the signature of 
the above admixture of the vortex free state. The dependence arises from the interference term. But the contribution 
of the latter becomes significant only after some time of duration of the process of free expansion. Taking the estimates 
of the root mean squared values of r±/R K and r\j ' R 2 K from the written spatial distribution function, one can find that 
the axial asymmetric contribution oc sin <fr becomes to dominate after the time of flight 

mRl I k 

Using the parameter sets Eqs. (||), [(§)], and the results of numerical solutions of Eq. @ together with Eq. JTsj), 
one can find from Eq. ( |i~4] ) that the interference term becomes to dominate after the free expansion time exceeding 
t = 24 ms (50 ms), respectively. Since the time of free expansion reported in Refs. [Q,[)| is 27 ms, it is clear that the 
proposed feature of the final wave function that could emerge after the stirring has not enough time to develop in 
the experiments . An additional testable feature of the proposed wave function of the final state is the flattening 
of the central vortex dip in the density distribution due to the increasing relative contribution of the zero angular 
momentum state oc 1 — k as time of free expansion of the BEC cloud is increasing. It would be interesting to enlarge 
(if possible) the time of ballistic expansion of the BEC clouds after the stirring to see if the dependence of their 
spatial distribution will acquire the angular dependence oc sin <f>, and to study the relative weights of the components 
of condensed atoms with different angular momentum at different times of their ballistic expansion. 
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FIG. 1. The dependence of the critical angular velocity of formation of the quantum state of BEC with one vortex on the 
portion k of atoms in the state with the unit angular momentum per particle. The curves labeled by A and B correspond to 
the trap parameters Eq. ([]) and (||), respectively. 
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